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Abstract. The exact dynamics of A'^ two-level atoms coupled to a common electromagnetic bath and 
closely located inside a lossy cavity is reported. Stationary radiation trapping eff'ects are found and very 
transparently interpreted in the context of our approach. We prove that initially injecting one excitation 
only in the A'' atoms-cavity system, loss mechanisms asymptotically drive the matter sample toward a 
long-lived collective subradiant Dicke state. The role played by the closeness of the A' atoms with respect 
to such a cooperative behavior is brought to light and carefully discussed. 

PACS. 03.65.Yz Decoherence, open systems, quantum statistical methods - 03.67.Mn Entanglement pro- 
duction and manipulation - 42. 50. Fx Cooperative phenomena in quantum optical systems 

1 Introduction ing entangled states of bipartite or multipartite systems. 

To produce and to be able to modify at will the degree 

It is well known that entangled states of two or more par- of entanglement stored in a system is indeed a desirable 

tides give rise to quantum phenomena that cannot be ex- target to better capture fundamental aspects of the quan- 

plained in classical terms. The concept of entanglement tum world. Over the last decade, moreover, it has been 

was indeed early recognized as the characteristic trait of recognized that the peculiar properties of the entangled 

the quantum theory itself. For this reason much interest states, both pure and mixed, can be usefully exploited as 

has been devoted by many physicists, both theoreticians an effective resource in the context of quantum informa- 

and experimentalists, toward the possibility of general- lion and computation processing [T1I21I51II) . Such a research 
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that quantum computers are, at least in principle, able to erning the dynamics of the system under scrutiny, suppos- 
solve very hard computational problems more efficiently ing that only one excitation has been initially injected in 
than classical logic-based ones. The realization of quan- it, we show that the system of the N two-level atoms may 
tum computation protocols suffers anyway the difficulty be guided, with appreciable probability, toward a nontriv- 
of isolating a quantum mechanical system from its cnvi- ial stationary condition described as a Dicke state having 
ronment. Very recently, however, nearly decoherence-free the form | S, —S) with 5* = ^^"^ , S being the total Pauli 
quantum gates have been proposed by exploiting, rather spin operator of the atomic sample. In addition, exploit- 
than countering, the same dissipation mechanisms [^]|(>]|7| ing the knowledge of the exact temporal evolution of the 
I51imil0lllllll2| . The main requirement to achieve this goal matter-cavity reduced density matrix, we propose an an- 
is the existence of a decoherence-free subspace for the sys- alytical route to follow up some physically transparent 
tem under consideration. In the same spirit it has been aspects characterizing the entanglement building up pro- 
recognized that transient entanglement between distant cess in the passage from a chosen totally uncorrelated ini- 
atoms can be induced by atomic spontaneous decay |13| tial situation to the manifestly entangled asymptotic one. 
or by cavity losses ^1] . In ref. ^l] it has also been demon- The treatment followed in our paper enables to catch the 
strated that asymptotic entangled states of two closely physical origin of the stationary collective Dicke behavior 
separated two-level atoms in free space can be created as of the system. In addition it has the virtue to provide a 
conseguence of the spontaneous emission process. transparent way to understand the key role played by the 

T ■ i ii 1 1-11 loss mechanisms and by the closeness of the atoms in the 

In this paper we present a new path along which loss 

, . ^ ^ ^. , . , . 11, • T^-i phenomena brought to light. The paper is organized as 

mechanisms act constructively inducing a collective Dicke o o t- t- a 

1 , . . li- i ITT n Pi n ^1 1 1, ■ , follows. The next section is devoted to an accurate pre- 

behavior m a multiatom sample. In reis. |15U1d| multistep '^^ 

1 , , irT-.-ij-i r A -L scutation of our physical model and to the formulation of 

schemes to generate a set of Dicke states of multi yl-type 

^1 , , , i 1 T ii 1 ii the relative master equation for the matter-cavity reduced 

three level atoms are reported. In these procedures the 

, • i • ii -1 -i-i r r 11 ■ i- density operator. An appropriate unitary change of this 

key point is the possibility of successfully incorporating j f t-t- t- j o 

,1 r -i 1 ■ ii ii 1 i- operator variable provides, in section 3, the mathematical 

the presence ol cavity losses m the theory, neglecting on t- i j 

^1 ^ , . ^ . . key tool for exactly solving a Cauchy problem in the one- 

the contrary atomic spontaneous emission. jo 

In order to reach our scope we consider a material sys- 
tem of N identical two-level atoms closely placed inside 
a resonant bad cavity taking also into account, from the 
very beginning, the coupling between each atom and the 
external world. Exactly solving the master equation gov- 



excitation subspace of N atoms-resonator Hilbert space. 
The two successive sections contain the main results of this 
paper. The entanglement formation process is addressed 
in section 4 studying the time evolution of the Wootters 
concurrence | 17II18| relative to a generic pair of two-level 
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atoms. Section 5, in turn, brings to light the occurrence physical system at i in a state with a well defined num- 
of stationary collective Dicke subradiant behaviour of our ber of excitations A^e, its dynamics is confined in a finite- 
matter subsystem. The last section contains some final dimensional Hilbert subspace singled out by this eigen- 
remarks as well as a discussion on the experimental im- value of N. In a realistic situation, however, the system 
plementation of the physical conditions assumed in the we are considering is subjected to two important sources 
paper. of decoherence. The first one is undoubtedly related to the 

fact that photons can leak out through the cavity mirrors 
due to the coupling of the resonator mode to the free radi- 
ation field outside the cavity. Moreover the atoms present 
inside the resonator can spontaneously emit photons into 
non-cavity field modes. The microscopic hamiltonian tak- 
ing into account these loss mechanisms may be written in 
the form ^ 



2 The physical system and its master 
equation 



As previously said, our system consists of N identical two- 
level atoms closely located within a single-mode cavity. 
Indicate the atomic frequency transition and the cavity 
mode frequency by ujq and u respectively and suppose H = H ac + Hn + Har + Hen. (2) 



where 



ojQ ^ uj. Assume, in addition, that all the conditions under 
which the interaction between each atom and the cavity .^-^ , , 

Hr = ?i2^t^k,A[Ck.ACk,A + Ck;^Ck,A] (3) 

field is well described by a Jaynes Cummings (JC) model, k,A 

is the hamiltonian relative to the environment, 

are satisfied T^. Thus, the unitary time evolution of the 

system we are considering is governed by the following Har = ^ [5k,A^k,Af + h.c] (4) 



hamiltonian: 



(0 ^(0 
k,A,i 

describes the interaction between the atomic sample and 



Hac - Wa + ;i^^f]a«+n^[£«a4^+/i.c.] (1) ^^^^ bath and, finally. 



2 

i=l 1=1 



HcR = ^[sk,ACk,Aa^ + h.c] (5) 



In eq. ^ a and denote the single-mode cavity field ^ ^ 
annihilation and creation operators respectively, whereas represents the coupling between the environment and the 
<j'i\ cTj!^ (i = I, ...N) are the Pauh operators of the «-th cavity field. In eq. Q we have assumed, as usual that 
atom. The coupling constant between the i — th atom and the two subsystems, the N atoms and the single-mode 
the cavity is denoted by e^*^. cavity, see two different reservoirs. In eq. Q- (0) the bo- 
It is easy to demonstrate that the excitation number son operators relative to the atomic bath are denoted by 
operator iV defined as TV = a^a+i ^^j^ cri*''-!-^ is a con- {ck,A,cj^>,} whereas Ck,A,Ck>, are the (k. A) mode anni- 
stant of motion being [N ,Hac] — 0. Thus, preparing the hilation and creation operators respectively of the cavity 
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environment. Moreover, the coupling constants {sk,A} are where 

phenomenological parameters whereas -ffis — ^ij'^t'^J (^) 



,sin(^r„) .,cos(^ri,),^ 
stems from a dipole atom-field coupling [21] ■ In eq. ||6ll (^Qrij) [ujQrijj'^ 

6kA represents the polarization vector of the atomic ther- CfPAc = k{2apAcct^ ~ ct^ctpAc ~ PAco:^o:) (10) 



mal bath (kA) mode of frequency uj^, V its effective vol- 
ume, d is the electric dipole matrix element between the 
two atomic levels and, finally, is the position of the i-th 



CaPac = (11) 



N 



ru(2iJL' PAC(^+ - (J+ crl'pAC - PACCryal') + 



i=l 

N 



atom. Indicate now by d and fi,- two unit vectors along the p /'o^(*) „ ^(j) „ „(»)„0')\ 



We point out that both the two baths entering in our 



atomic transition dipole moment and the atomic distance 
Tij = r,; — Vj, respectively. Our microscopic hamiltonian 
model (|2Jl does not take into account dipole-dipole static 
interactions among the N atoms making in this way H in- 
variant under the exchange of two arbitrary atoms when 
Tij << ^ for any i and j. The two hypotheses imply- 
ing such a permutational symmetry property are, in gen- 

eral, contradictory ^E^. However, notwithstanding the Moreover in eq. ^ the couphng constants 
conceptual difficulties connected with these assumptions, 47r?ia;'^|dP 
they are commonly adopted |l!'{lll4ll24ll25j for the sake of 
simplicity. ~ ~ * 7^ J (14) 



model are supposed in thermal states at T = and that, 
when Tij << f2ij tends toward a limiting value [2^ 
independent from i and j |28| . 

The decay rate k appearing in eq. pU|l is given by 

/c = ^|skA|'5(wk-w) (12) 

kA 



r.^r.^i^ (13) 



Following standard procedures based on the Rotating- related to the spontaneous emission loss channel, define 

Wave and Born-Markov approximations fMF7\ . it is pos- tl^c spectral correlation tensor T 

sible to prove that the reduced density operator pAc rel- In eq. ^ the function f,j is defined as follows: 

ative to the bipartite system composed by the N atoms f _ ^rn - ^2l sin(^rij) 

I LOorij 

subsystem and the single-mode cavity, evolves nonunitar- cosf— r ) sinf— r- ~) 



ily in time in accordance with the following Lindbland 



[1 - 3(d . r,,r]\c^ y^ c^-^^]} 



master equation |2H1' important to underline that the last term appearing 

in the right hand side of equation is a direct con- 

PAC = -^[Hac + HlStPac] + C,fPAC + Cap AC (7) sequence of the fact that we have considered, from the 
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Hac= U^HacU^ (19) 

N 

hioa^a + fi-^ ^ cr^') + fee// [a(T+ -* + h.c] 



very beginning, a common bath for the N atoms. As we for i — I, N so that 
shall see, this term is responsible for cooperative effects 
among the N atoms leading to the possibility of generating 

asymptotic entangled states of the atomic sample, immune ^ i=i 

from decoherence. We wish to stress that it is the nearness ^i*^ ^-ff = Transforming in the N, = 0,1 excita- 

of the atoms inside the cavity that imposes the considera- ^ion subspace the operator variable pAc into the new one 

tion of a common bath. If, on the other hand, the distance ~PAC = W PAcU and taking into account that 
among the atoms became large enough {nj > ^), these jj^^ = U'^HlsU = HlNct^I'' <t^1\ (20) 

cooperative effects, as deducible from eq. (|15|) . would dis- 



being the common limiting value of Qij when tends 
to zero, it is not difficult to convince oneself that 



appear so that the dynamics of the system could be equiv- 
alently obtained considering N different reservoirs, one for 

each two-level atom. In such a situation the system would Pac — "'^['^-^'^ ^ H^s, Pac] (21) 

evolve toward its vacuum state with no excitation. + k{2apAcCi^ — Q.^ <xpAC — PAcct^ o;) 

, (1), (1) (1) (1). . (1) (1). 
+ Nr{2al'pAc<7+ - cr+ crl'pAc - PAC(^+ ) 

3 One-excitation exact dynamics ■ ■ r m rm rrm mrrr ^i 

m view of eq. lfIo|) - If20|) . Comparmg eq. (|2Tjl with 

In what follows we assume that the atoms within the cav- ^q. shows that in the new representation the corre- 

ity are located at a distance smaller than the wavelength spondent spectral correlation tensor f is in diagonal form, 

of the cavity mode, thus legitimating the henceforth done moreover being Pu = Fdi,,. The physical meaning of this 

position £(*) = e and g[:^^ = gk,x for any i. Under this hy- Pecuhar property is that the atomic subsystem in the 

pothesis we solve eq. © exploiting the unitary operator transformed representation looses its energy only through 

U defined as interaction of the first atom with both the cavity mode 



N 



U = Y[U, (16) 

where 

g5,(4<TL-<Ti^;) i = 2,...,N (17) 

with 6i = — arctan( ) and [C/, TV] = 0. It is easy to 
demonstrate that, if no more than one excitation is ini- 
tially stored in the atom-cavity physical system 



and the environment. Such a behaviour stems from the 
fact that, in view of eq. (|19|l . the other — 1 atoms freely 
evolve being decoupled either from the cavity field and 
from the electromagnetic modes of the thermal bath. It 
is of relevance to underline that the form assumed by the 
terms associated to the nonunitary evolution, appearing 
in eq. 1)21(1 . directly reflects the main role played by the 



C/V^-'C/ = \ - — -cr^^ — . / ——^ — -a^^ H ^(T^P . closeness of the atoms in our model. It is indeed just this 

V i + y 1(1 - 1) + jn + 



l{l-l) 

(18) feature which leads, in the transformed representation, to 
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the existence of A^ — 1 collective atoms immune from spon- 
taneous emission losses and, at the same time, decoupled 
from the cavity mode. Thus, to locate the atomic sample 
within a linear dimension much shorter than the wave- 
length of the cavity mode, introduces an essential per- 
mutational atomic symmetry which is at the origin of a 
collective replay of the iV atoms such that, even in pres- 
ence of both the proposed dissipation routes, the matter 
subsystem may stationarily trap the initial energy. 

Bearing in mind that [Hac,N] = and iV] — 0, 
it is immediate to convince oneself that, if only one ex- 
citation is initially injected into the atomic subsystem, 
whereas the cavity is prepared in its vacuum state, at 
a generic time instant t, the density operator paCj can 
have not vanishing matrix elements only in the Hilbert 
subspace generated by the following ordered set of + 2 
state vectors: 

|0)|-)lh)2...|-)A.^|/3l) 

|0)|-)i...|+>^...|-)jv = |/3h+i) h^l,..,N 
|l)h)i...h)^^ 1/3^+2) 

where \p) {p = 0, 1) is a number state of the cavity mode 
and (|— denotes the excited (ground) state of the 
h-th collective atom {h = 1,...,N). Eq. (|21(l can thus 
be easily converted into a system of coupled differential 
equations involving the density matrix elements jDhj = 
{(3h\p\Pj) with h,j = 1,...N + 2. At this point let's ob- 
serve that from an experimental point of view it seems 
reasonable to think that the excitation given at t = to 



the matter sample can be captured by i-th the atom or 
by j-th with the same probability. In other words our ini- 
tial condition must be reasonably represented as statistical 
mixture of states \Ph), with ft, > 2, of the form 



N+l 



(22) 



h=2 



Exploiting eq. (fT5|) it is possible to verify that 



(23) 



After lengthly and tedious calculations and taking into 
account eqs. H22|l - H23(l . we have exactly determined the 
time evolution of each p^j {t) h,j = 1, ..,N + 2 finding: 

PAc{t) = 
I 



Pii(t) 

P22(t) 

P33(t) 

m 



P2,N + 2{t) 



(24) 



... PN + l,N + l(t) 
y P2.N + 2W ■ ■■■ PJV + 2,JV + 2(t) J 

where p,j(t) = p,j(0) {3 < i,j < N + 1), pi,i(t) = 1 
ES^Pi.»(0 and 

-{k+Nr)t 



P22{t) = 



[(a^ + 6^ + |Z\|^)cosh(&i)-)- 



2N{a? + fe2) 
(a^ + _ |zi|2)cos(at) - 

2(6(0)0 — 1^) — aj4_) sin(at) + 

2(a(a)o - oj) + M_) sinh(fet)] 



(25) 



P2,N+2{t) 



^effe 



-{k+Nr)t 



PN+2,N+2{t) 



-|-Z\(cosh(M) - cos(at))] 
2^2 g-(fc+7vr)t 



[Q{i sin(ai) + sinh(6i)) 



(26) 

cosh(6i) - cos(ai)] (27) 
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with of how to extend to a generic A^-partite physical sys- 

-l tern definition and measure of entanglement built for bi- 

a ={-[{^,-uf + Aelff-A^ (28) 

partite systems, constitutes a topical challenge involving 

+ \[{{^,^u;f + Aelff-A,f+A{Q,-ufAlY}K ^^^^^^^ 
^ many researchers |,i()U,ilU,i/U.i,iU,i4U,SnU,inj . is is well under- 

stood from first principles that when many subsystems of 
^ 2 ^^"^ ^ '^^^ ^-^ ^ ^ a multipartite system individually entangle a prefixed one, 

\ 2 2 2 2 2 -1 

+ 2 [(('^0 — w) + 4egyj — A-) + 4(a)o — lo) A_'\ }^ , ^^^^ entanglement degree within each pair, anyhow mea- 
sured, is subjected to quantitative restrictions [HOIESIEI 
I35| . This, for instance, implies that two maximally entan- 
gled subsystems of a multipartite system are necessarily 
disentangled from any other constituent units of the to- 
tal system in an arbitrarily given pure or not state. Of 
course, whatever the multipartite entanglement definition 
is adopted, its occurrence is conceptually compatible with 
a complete lack of partial entanglement of a given order 
for example binary. On the contrary the existence of en- 
tanglement between two specific subsystems has to be con- 
sidered as a clear symptom of entanglement in the mul- 
tipartite system. Following this line of reasoning we here 
therefore propose to study the time evolution of the en- 
tanglement within all the possible binary subsystems that 



and Q = a + ib = ^[{ujq-uj)^ -A'i+4:slj^+iA^{Cjo~u;)], 
A^ = k - NT, A ^ ujQ - u! + lA^ and Cjq ^ ojq + ^J7. 
We wish to emphasize that, on the basis of the block di- 
agonal form exhibited by paCj a generic time instant 
the transformed matter-radiation system is in a sta- 
tistical mixture of its vacuum density matrix and of an 
one-excitation appropriate density matrix describing with 
certainty the storage of the initial energy. Eqs. H25() - (|27() . 
giving the explicit form of the time evolution of the com- 
bined physical system, allow the exact evaluation of the 
mean value of any physical observable of interest and, for 
instance, to follow the entanglement formation or the pro- 
gressive raising up of decoherence effect in the matter- 
cavity subsystem. 



4 Entanglement building up 

The circumstance that we succeed in finding the explicit 
time dependence of the solution of the master equation 
H21() . provides a lucky and intriguing occasion to analyze 
in detail at least some aspects of how entanglement is 
getting established in our exemplary enough multipartite 
system. We wish indeed to point out that the question 



is within each of the ^'^^^^^ pairs of individual parts ex- 
tractable from the iV-partite set under scrutiny. To this 
end we choose to evaluate the Wootter's concurrence [171 
I18| Cij (t) to characterize quantitatively the formation of 
entanglement within the —th pair of two-level atoms 
of our matter sample. Since the dynamical problem ex- 
actly solved in the previous section is invariant under the 
exchange of two arbitrary atoms, then one may guess and 
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indeed easily prove, that the reduced density matrix +PN+2.N+2{t) (35) 

N -2 N -1 

= TrA,^{UpAU^ (30) b = c^ ——~p^^{t) + ——~p^^^{t) (36) 

r .1 /■ •\ • • ^. + n ■ A A d = e= ^{p22{t)- PN,N{t)) (37) 

oi the (I, J) — tn pan is structurally independent irom A* 

the indices of two prefixed atoms meaning that the sub- It is now easy to construct and diagonahze R^'^~'^'^\t) 

stitution of i and j with i' and / respectively, exactly finally getting for any pair with i < j 

maps into p^'^'l The symbol Tta,, means to trace 2\p22{t) - pn+i N+i{t)\ , . 

over the atomic variables excluding the pair (i, j) whereas 2^i=2 PiH'-J 

~ _ 5 ?, \ The c— function Cij{t) is displayed in figures (l)-(3) in 

Pa — ' cavity \P AC S ■ 

The concurrence Q,(i) is defined as correspondence to A^ = 2, 10 and 100 respectively us- 

, — ■ — , — ■ — I — ■ — I — ■ — ing reasonable values for the involved parameters 1^1 ■ It 

Q,(t) = max(0, VaI^^'- Va^'^'^- VAf VAf'^ (31) 

represents the conditional concurrence characterizing the 

where aI'"''' (g = 1, 4) are the decreasing-ordered eigen- , , , ,. r u- -t-j- -t i ^ a t-i, 

•J ' ' ' o t) temporal evolution ot bipartite entanglement under the 

values of the matrix , -ti, • ii, ^ ii, i t, -t • v n • • 4- j ■ j-u 

hypothesis that the only photon, initially injected m the 



_, - s , fact that Cij{t) is different from zero at any time instant 



i?^*-') p^*-?) • (32) system has not escaped because of loss mechanisms. The 

where the spin flipped matrix p^'-'^ is given by 

t > whatever the pair is, undoubtedly reflects the ex- 

=(T(^)®a(j')(pfe))VW®4^^ (33) 

istence of a process giving rise to the entanglement for- 
(pf*-?))* being the conjugate matrix of [TTItTH]. mation inside the A^-partite system. On the other hand 

In view of the invariance property of p'-»\ it is not dif- since each atom turns out to be entangled with aU the 
ficult to persuade oneself that the eigenvalues of i?(*-')(t) others, the concurrence of each couple of atoms is mono- 
are pair-independent, which, as a consequence, implies tonically decreasing with A^. As already mentioned at the 
C'ij(i) — CN-i,N{t) for any i,j — 1,2, ...,N {i < j). Let's beginning of this section, this behaviour reflects nothing 
then start by extracting the expression of p'^^-^'^^t). but the expected reduction of atom-atom entanglement 
Tracing in accordance with eq. (jSUJ) yields due to the increase in the number of possible entangled 

/,\ I \ I \ / I / I /o/i\couples. Albeit Caii) tends to vanish when the number 

PN-l,N{t) = a\ — ) N-l\ — ) N N-i'\ — \n'\~\ io4j ^ "-J^ ' 

, , , , , , , , , , A^ of atoms goes to infinity, we find the remarkable result 

+ &|+)Ar-l|-)Ar W-l(+|Ar(-| +c|-)Ar_i|+)Ar Ar_i(-|Ar(+| 

that the total binary concurrence 

+ d\^)N-l\~)N W-l(-U( + | +e|-)Ar_i|-|-)jv Ar^l( + |jv(-|, 

CBT{t)^Y.^^Ai) = 

where 

2 2 - - 1)IP22(^) - PN+l,N+l{t)\ 

a=il--)p22{t) + iN-3 + -)p^Mt)+ - Elfpuit) ' ^ ^ 
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1 2 3 4 5 



Fig. 1. Conditional binary concurrence Cij{t) in correspon- 
dence to iV = 3, £ = KfHz, k = lO'^Hz, F = lO^Hz and 

UJ = Coo= 10"H2 



exhibits an oscillatory time behaviour with a decreasing 
amplitude around a time-dependent mean value mono- 
tonically tending toward the stationary value 1, whatever 
N is (see figures (4)-(6)). This dynamical property is of 
relevance because examples of multipartite system states 
manifestly entangled, for which Cbt = may be provided 
|34j . Stated another way, eq. (|39|l tell us that the not 
vanishing contribution of the total binary entanglement 
to the formation of entanglement within our multipartite 
dynamical problem does not scale with N , reflecting that 
the A^— decrease of each Cij{t) is well compensated by the 
quadratic A'^— increase of the number of pairs. We thus 
claim that the behaviour of CsTit), and in particular its 
asymptotic tendency, provides, in our situation, a peculiar 
feature helping to describe and understand some aspects 
of the entanglement formation in our multipartite system. 



Ci,(t) 




0.5 1 1.5 2 2.5 xlO ^sec 

Fig. 2. Conditional binary concurrence Cij(t) in correspon- 
dence to iV = 10, e = ICf'Hz, k = W^Hz,r = lO^Hz and 
tj = tJo = W^'^Hz 

Ci,(t)r 



.0002 
.00015 

. 0001 
0.00005 




0.1 0.2 0.3 0.4 xlO'^sec 

Fig. 3. Conditional binary concurrence Cijit) in correspon- 
dence to iV = 10, e = lO^Hz, k = W^Hz,r = lO^Hz and 
lo^lSq = 10^'^ Hz 



1 




0.1 0.2 0.3 0.4 0.5 ms 



Fig. 4. Total binary concurrence CBrit) in correspondence to 
iV = 3, e = IQ^Hz, k = 10^Hz,r = lO^Hz and w = wo = 
lO^^Hz 
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xlO sec 

Fig. 5. Total binary concurrence CsTit) in correspondence to 
iV = 10, £ = lO^Hz, k = W^Hz,r = lO^Hz and lu ^ cJo ^ 

1 

. 8 
. 6 
0.4 
0.2 




. 1 



. 2 



. 3 



0.4 



. 10 sec 



Fig. 6. Total binary concurrence CBxit) in correspondence to 
iV = 100, e = IQ^Hz, k = 10*Hz,r = lO^Hz and = oJo = 
10"J/2 

5 The asymptotic form of pac 

It is of particular relevance that for 



t> {k + NT)-^ =TAC 



(40) 



t ^ TAG, the eq. H24|l assumes the following diagonal 
asymptotic form 



PAc{t > TAC) 



^^0 




0^ 












^0 




(41) 



Starting from eq. H41() and coming back to the old rep- 
resentation, it is possible to give at any time instant t the 
exact solution pAc for the reduced density matrix of the 
system under scrutiny. Taking into account that the uni- 
tary operator U is time independent, in view of eq. H41|l 



PAC = UpAc{t > TAc)U^ 



(42) 



is time independent too. In fact we find that for t 3> tac, 
the reduced density matrix can be written in the compact 
form 

N+l 



PAC = + ^ E I^T ^)(^' 

h=2 



(43) 



1 



l/3i)(/3i 



N - 1 



N 



-Pn 



where 



N+l 



A) 



(44) 



(h) 



the correspondent asymptotic form assumed by pAc is It is worth noting that each normalized state jc^c^ 

time independent and such that the probability of find- given by eq. (|44() . defines a particular subradiant Dicke 

ing energy in the effective JC subsystem exactly vanishes, state. It is indeed possible to prove that, whatever N and 

Taking into account the easily demonstrable inequality h are, the states l^s^''') are common eigenstates of Sz and 



b < Tjjj, it is immediate to convince oneself that for S- pertaining to eingevalues 



JV-2 



and respectively. 
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This property is a sufficient condition to claim that l^'^ ^> ^ Conclusive Remarks 

is eigenstate of too, having the form \S, —S) with S — 

Summing up, in this paper we have exactly solved the 

^^Y^. This means that 

dynamics of N identical atoms resonantly interacting with 

S^PrnixS^ = (5(5 + l))^pmix (45) a single mode cavity, taking into account from the very 

beginning the presence of both the resonator losses and the 

S zPmixS z — 5 Pmix (46) 

atomic spontaneous emission. We have moreover supposed 

and that each IV't''^) defines an example of subradiand or that only one excitation is initially injected into the system 

trapped state [57113811391140; . Thus the result expressed by of interest and that the atoms are located in such a way 

eq. 1031) suggests that a statistical mixture of stationary to experience the same cavity field. 

subradiant Dicke states of the atomic sample, having well From a mathematical point of view the novelty of our 
defined values of and 5^, can be generated, at least results is the presentation of an exact way to solve the 
in principle, putting outside the cavity single photon de- master equation eq. Q of the system, based on the uni- 
tectors allowing us to continuously monitor the decay of tary transformation accomplished by the operator U given 
the system through the two possible channels (atomic and by eq. (dnj. In the new representation associated to such 
cavity dissipation) Eq. (gSJ, indeed, clearly shows a specific [/, the differential equations governing the tem- 
that, reading out the detectors state at f > tac, if no poral behavior of the density matrix elements p/jj , become 
photons has been emitted, then, as a consequence of this indeed much simpler to solve if compared with the ones 
measurement outcome, our system is projected into the ruhng the temporal behavior of phj. This circumstance 
state jfJ2h=2 \''l'T ^)i''pT ^ \ - stems from the fact that in the new representation only 
Stated another way, successful measurements, performed one atom is at the same time coupled both to the cavity 
at large enough time instants t, generates an uncorrelated field and to the quantized electromagnetic modes of the 
state of the two atomic and cavity field subsystems, leav- thermal bath. The decoupling of the — 1 atoms is a di- 
ing the matter sample in the statistical mixture of Dicke rect consequence of the permutational symmetry proper- 
states satisfying eqs. and (^BJ. On the basis of ties acquired by the matter subsystem under the assumed 
the analysis reported in the previous section it is possible point-like model condition. 

to state that such a statistical mixture is entangled. From a physical point of view the new results reported 

Let's finally observe that the probability P that at in this paper have the merit of providing the key for trans- 

t ~ t the excitation is still contained in the atomic system parently interpreting the origin of the asymptotic occur- 

increases with the number N of atoms, being P = 1 — jj, rence of collective subradiant Dicke behavior of the matter 

as immediately deducible from eq. (|43(l . subsystem. The analysis and the discussion presented in 
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section 4, highlighting some features of the entanglement of eq. (|40|l . the value of tag correspondent to fc = W^Hz 
formation process, legitimate the claim that the asymp- and F ~ lO'^Hz |S] becomes much less than 10^^ what- 
totic condition toward which our physical system is guided ever N is. The experimental implementation of the spe- 
by the loss mechanisms exhibits entanglement. It is of rele- cific conditions envisaged in our paper thus require the 
vance to notice that the form of the stationary conditional ability of locating for a time of the order of 10~^sec, N 
state Pmix appearing in eq. H43|l is independent from both atoms in an enough small region within an optical res- 
r and fc, while the exponential tendency toward the sta- onator. In particular the distance between two arbitrarily 
tionary condition is governed by the rate k + NT. This chosen atoms of our matter sample has to be much less 
means that just the presence of only one loss channel is suf- than 500nm. The more and more growing technological 
ficient to address the same asymptotic radiation trapping successes registered in the last few years in the confine- 
condition, even if the transient duration is characterized ment of individual atoms j45ll46ll47ll48| or clouds of iden- 
by a different time decay constant. We wish in addition to tical atoms |45ll4(j|l49| with high spatial resolution, suggest 
emphasize that if the microscopic model neglects sponta- that implementing our conditions in the near future is in 
neous atomic decay, cooperative effects occur even if the the grasp of the experimentalists. To enforce our claim it 
atomic sample is spatially dispersed j41ll42ll4l-{ll44j . On the is appropriate and relevant to quote the paper of Ozeri 
contrary, the closeness of the N atoms is a necessary re- et al wherein the authors experimentally demonstrate 
quirement when the complete Hamiltonian model Q is the possibility of confining in a blue-detuned optical trap 
used (regardless of how bad the cavity is) in order that a a sample of 10^ Rubidium atoms in a region of 0.1-^30nm 
robust entanglement may be conditionally reached in the for O.Ssec at a temperature of 24/i_?C. 
matter subsystem. If indeed the distance among the N 
atoms is larger than the radiation wavelength, collective Rgfgrgficgs 
behavior stemming from the interaction of each atom with 

1 ,1 J. .,1 ,1 ,1 , ,1 1. A.Beiee, Hueo Cable, P.L.Knieht, in Proceedins of SPIE, 

a common bath, disappear with the consequence that the & i & i & i a ; 



probability of finding in the system the initial energy goes 
toward zero with time. We thus may state that the key for 
trapping the energy in the atomic sample, inducing a sta- 
tionary collective Dicke behavior, is the closeness among 
the atoms. We wish to conclude presenting some remarks 
concerning the experimental relevance of the problem dis- 
cussed in this paper. We begin by observing that, in view 
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